Using a Barycenter Compact Model for a Circuit Network

ABSTRACT

Any primitive cells or blocks can be represented physically by a Barycenter compact model (or Barycenter model), and any black box model can also be physically represented by a Barycenter compact model physically. A hierarchical boundary condition between blocks is formulated by the Barycenter compact model. Hierarchical boundary condition problems between blocks can be limited within two levels only if using the Barycenter compact model.

CROSS-REFERENCE TO RELATED APPLICATIONS

This patent application is a continuation of U.S. patent applicationSer. No. 14/248,269, filed Apr. 4, 2014, issued as U.S. Pat. No.9,111,058 on Aug. 18, 2015, which is a continuation of U.S. patentapplication Ser. No. 13/178,288, filed Jul. 7, 2011, issued as U.S. Pat.No. 8,694,302 on Apr. 8, 2014, which is a continuation-in-partapplication of U.S. patent application Ser. No. 13/159,384, filed Jun.13, 2011, which claims the benefit of U.S. provisional patentapplication 61/354,186 and 61/354,189, both filed Jun. 11, 2010, and isa continuation-in-part of U.S. patent application Ser. No. 12/915,362,filed Oct. 29, 2010, which is a continuation of U.S. patent applicationSer. No. 11/421,206, filed May 31, 2006, issued as U.S. Pat. No.7,827,016 on Nov. 2, 2010, and Ser. No. 11/421,212, filed May 31, 2006.These applications are incorporated by reference along with all otherreferences cited in this application.

BACKGROUND OF THE INVENTION

This present invention relates to the field of electronic designautomation for electronic circuits, and more specifically, to systemsand techniques to solve a network using a Barycenter compact model and ahierarchical scheduler.

The age of information and electronic commerce has been made possible bythe development of electronic circuits and their miniaturization throughintegrated circuit technology. Integrated circuits are sometimesreferred to as “chips.” Some types of integrated circuits includedigital signal processors (DSPs), amplifiers, dynamic random accessmemories (DRAMs), static random access memories (SRAMs), erasableprogrammable read only memories (EPROMs), electrically erasableprogrammable read only memories (EEPROMs), Flash memories,microprocessors, application specific integrated circuits (ASICs), andprogrammable logic

Integrated circuits have been widely adopted and are used in manyproducts in the areas of computers and other programmed machines,consumer electronics, telecommunications and networking equipment,wireless network and communications, industrial automation, and medicalinstruments, just to name a few. Electronic circuits and integratedcircuits are the foundation of the Internet and other on-linetechnologies including the World Wide Web (WWW).

There is a continuing demand for electronic products that are easier touse, more accessible to greater numbers of users, provide more features,and generally address the needs of consumers and customers. Integratedcircuit technology continues to advance rapidly. With new advances intechnology, more of these needs are addressed. Furthermore, new advancesmay also bring about fundamental changes in technology that profoundlyimpact and greatly enhance the products of the future.

To meet the challenges of building more complex and higher performanceintegrated circuits, software tools are used. These tools are in an areacommonly referred to as computer aided design (CAD), computer aidedengineering (CAE), or electronic design automation (EDA). There is aconstant need to improve these electronic automatic tools in order toaddress the desire for higher integration and greater complexity, andbetter performance in integrated circuits.

Large modern day integrated circuits have millions of devices includinggates and transistors and are very complex. As process technologyimproves, more and more devices may be fabricated on a single integratedcircuit, so integrated circuits will continue to become even morecomplex with time. In the past, many parasitic effects may not have beenconsidered because they were less significant or insignificant comparedto other factors.

As lithography and miniaturization techniques advance, on-chip devicesand line widths become smaller, frequencies increase. As a consequence,many more impedances such as parasitic resistances, inductances, andcapacitances and parasitic effects need to be considered. If theseparasitics and effects are not taken into account, poor simulationresults will result, and possible the electronic circuits will not workas expected after the circuit is fabricated. As more and more parasiticand other effects are accounted for, the circuit networks to besimulated become much more complex. As complexity increases, simulatingthe network takes significantly more computing resources and computationtime.

More specifically, in nanometer, gigahertz, low power VLSI design,power, and signal integrity has become critical. To accurately analyzechip performance, it is desirable to consider the impact of powerfluctuation, and the capacitive, inductive, or even substrate couplingnoise with devices, or any combination of these. This analysis entailsconsidering a very large amount of elements, which results in a verylarge system matrix for circuit simulation. This is a lack of a circuitsimulation algorithm that can simultaneously resolve a large number oflinear or linear with nonlinear devices while maintaining bothefficiently and accuracy.

Some problems with the prior art are performance, capacity (millions ofelements), accuracy (iterative matrix solver may be divergent for largenetwork and hard for parallel processing and distributed computing), thesize of memory has limitation in the computer, multi-thread is limitedby memory size, distributed computing (does not share memory),diakoptics (tearing and reassembly, but no reassembly method forhierarchical design), and no efficient way to formulate the hierarchicalboundary condition and solve the problem.

Therefore, there is a need for tools for solving networks.

BRIEF SUMMARY OF THE INVENTION

A technique uses a Barycenter compact model and hierarchical schedulerto solve a large network in a computing environment. The computingenvironment can be homogeneous or heterogeneous such as multiple cores,grid, and networking together. A tool of the invention solves a circuitnetwork in a hierarchical, multicore, and distributed computingenvironment and obtains an exact solution.

In an implementation, memory usage is controlled by I/O slicing to fit alimitation of the computer memory system. In another implementation,primitive cells or blocks are represented physically by a Barycentercompact model (which may also be referred to as the Barycenter model). Ablack box model can also be physically represented by a Barycentercompact model physically.

In an implementation, a hierarchical boundary condition between blocksis formulated by the Barycenter compact model. Hierarchical boundarycondition problems between blocks can be limited within two levels ifusing the Barycenter compact model.

A method is provided for solving a boundary condition of a large networkby a direct method. A method is provided for solving a boundarycondition of a large network by an iterative method. A method isprovided for solving a boundary condition of a large network by both adirect and an iterative method.

In an implementation, a hierarchical scheduler is generated by ahierarchical netlist partition or user defined. In anotherimplementation, job dependence and job submittal in a multicore ordistributed computing environment is managed by a hierarchicalscheduler.

Some Concepts include primitive cell; block; branch, meaning impedanceor admittance, such as R (resistance), C (capacitance), L (inductance)or Mutual Impedance such as M, or K; I/O branch current; I/O branchvoltage; internal branch current; internal branch voltage; nodalvoltage; mesh current; mesh voltage; ideal current source; ideal voltagesource; full pins equivalent model; I/O equivalent model; black boxmodel; Barycenter compact model, a Barycenter is the center of mass oftwo or more bodies which are orbiting each other, and is the pointaround which both of them orbit; Barycenter compact model; multiplecores, Networking and Grid Computers; and hierarchical database.

Other objects, features, and advantages of the present invention willbecome apparent upon consideration of the following detailed descriptionand the accompanying drawings, in which like reference designationsrepresent like features throughout the figures.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows a system of the present invention for using a Barycentercompact model and hierarchical scheduler to solve a large network inhierarchical, multicore, and a distributed computing environment.

FIG. 2 shows a system block diagram of computer system 101 used toexecute software of the present invention.

FIG. 3 shows a black box model.

FIG. 4A shows a first Barycenter compact model block with an internalcenter point.

FIG. 4B shows a second Barycenter compact model block with a boundarynode as a center point.

FIG. 5A shows a black box model being replaced with a first Barycentercompact model with an internal center point.

FIG. 5B shows a black box model being replaced with a second Barycentercompact model with a boundary node as a center point.

FIG. 6 shows an implementation of interconnected blocks within a block.

FIG. 7 shows interconnected blocks with a Barycenter compact model.

FIG. 8 shows a block diagram of a Barycenter compact model.

FIG. 9 shows a graphical diagram of a system environment input.

FIG. 10 shows a flow for a bottom up calculation.

FIG. 11 shows a flow for a top down update.

FIG. 12 shows a flow of an output.

FIG. 13 shows a specific implementation of a hierarchical partition.

FIG. 14 shows a hierarchical tree scheduler.

FIG. 15 shows a hierarchical dependence scheduler.

FIG. 16 shows a large I/O slicing into smaller I/O blocks.

FIG. 17 shows a flow diagram of a static simulation.

FIG. 18 shows a flow diagram of a dynamic simulator.

FIG. 19 shows a flow diagram of a circuit simulator.

FIG. 20 shows an example of a branch of a circuit.

FIG. 21 shows a graph of a tree and links.

FIG. 22 shows local trees, global links, and local links of a graph.

DETAILED DESCRIPTION OF THE INVENTION

FIG. 1 shows a system of the present invention for using a Barycentercompact model and hierarchical scheduler to solve a large network inhierarchical, multicore, and a distributed computing environment. In anembodiment, the invention is software that executes on a computerworkstation system, such as shown in FIG. 1. FIG. 1 shows a computersystem 101 that includes a monitor 103, screen 105, cabinet 107,keyboard 109, and mouse 111. In some implementations, however, themodeling and simulating systems are operated headless, which means thesesystems will not a mouse, keyboard, display, and so forth. Mouse 111 mayhave one or more buttons such as mouse buttons 113. Cabinet 107 housesfamiliar computer components, some of which are not shown, such as aprocessor, memory, mass storage devices 117, and the like.

Mass storage devices 117 may include mass disk drives, floppy disks,magnetic disks, optical disks, magneto-optical disks, fixed disks, harddisks, CD-ROMs, recordable CDs, DVDs, recordable DVDs (e.g., DVD-R,DVD+R, DVD-RW, DVD+RW, HD-DVD, or Blu-ray Disc), flash and othernonvolatile solid-state storage (e.g., USB flash drive),battery-backed-up volatile memory, tape storage, reader, and othersimilar media, and combinations of these. A binary machine-executableversion of the software of the present invention may be stored or resideon mass storage devices 117. Furthermore, the source code of thesoftware of the present invention may also be stored or reside on massstorage devices 117 (e.g., magnetic disk, tape, CD-ROM, or DVD).

A computer-implemented version of the invention may be embodied using,or reside on, computer readable medium. A computer-readable medium mayinclude any medium that participates in providing instructions to one ormore processors for execution. Such a medium may take many formsincluding, but not limited to, nonvolatile, volatile, and transmissionmedia. Nonvolatile media includes, for example, flash memory or opticalor magnetic disks. Volatile media includes static or dynamic memory,such as cache memory or RAM. Transmission media includes coaxial cables,copper wire, fiber optic lines, and wires arranged in a bus.Transmission media can also take the form of electromagnetic, radiofrequency, acoustic, or light waves, such as those generated duringradio wave and infrared data communications. Storage and computing maybe via the Internet such as on the cloud or cloud computing.

For example, a binary, machine-executable version, of the software ofthe present invention may be stored or reside in RAM or cache memory, oron mass storage device 117. The source code of the software of thepresent invention may also be stored or reside on mass storage device117 (e.g., hard disk, magnetic disk, tape, or CD-ROM). As a furtherexample, code of the invention may be transmitted via wires, radiowaves, or through a network such as the Internet.

FIG. 2 shows a system block diagram of computer system 101 used toexecute software of the present invention. As in FIG. 1, computer system101 includes monitor 103, keyboard 109, and mass storage devices 117.Computer system 101 further includes subsystems such as centralprocessor 202, system memory 204, input/output (I/O) controller 206,display adapter 208, serial or universal serial bus (USB) port 212,network interface 218, and speaker 220. The invention may also be usedwith computer systems with additional or fewer subsystems. For example,a computer system could include more than one processor 202 (i.e., amultiprocessor system) or the system may include a cache memory.

The processor may be a dual core or multicore processor, where there aremultiple processor cores on a single integrated circuit. The system mayalso be part of a distributed computing environment. In a distributedcomputing environment, individual computing systems are connected to anetwork and are available to lend computing resources to another systemin the network as needed. The network may be an internal Ethernetnetwork, Internet, or other network. Some examples of distributedcomputer systems for solving problems over the Internet includeFolding@home, SETI@home, and the Great Internet Mersenne Prime Search(GIMPS).

Arrows such as 222 represent the system bus architecture of computersystem 101. However, these arrows are illustrative of anyinterconnection scheme serving to link the subsystems. For example,speaker 220 could be connected to the other subsystems through a port orhave an internal connection to central processor 202. Computer system101 shown in FIG. 1 is but an example of a computer system suitable foruse with the present invention. Other configurations of subsystemssuitable for use with the present invention will be readily apparent toone of ordinary skill in the art.

Computer software products may be written in any of various suitableprogramming languages, such as C, C++, C#, Pascal, Fortran, Perl, MatLab(from MathWorks, Inc.), SAS, SPSS, Java, JavaScript, and AJAX. Thecomputer software product may be an independent application with datainput and data display modules. Alternatively, the computer softwareproducts may be classes that may be instantiated as distributed objects.The computer software products may also be component software such asJava Beans (from Oracle) or Enterprise Java Beans (EJB from Oracle).

An operating system for the system may be one of the Microsoft Windows®family of operating systems (e.g., Windows 95, 98, Me, Windows NT,Windows 2000, Windows XP, Windows XP x64 Edition, Windows Vista, Windows7, Windows 8, Windows CE, Windows Mobile), Linux, HP-UX, UNIX, Sun OS,Solaris, Mac OS X, Apple iOS, Android, Alpha OS, AIX, IRIX32, or IRIX64,or combinations of these. Other operating systems may be used. Eachcomputer in a distributed computing environment may use a differentoperating system.

Furthermore, the computer may be connected to a network and mayinterface to other computers using this network. For example, eachcomputer in the network may perform part of the task of the many seriesof circuit simulation steps in parallel. Furthermore, the network may bean intranet, internet, or the Internet, among others. The network may bea wired network (e.g., using copper), telephone network, packet network,an optical network (e.g., using optical fiber), or a wireless network,or any combination thereof. For example, data and other information maybe passed between the computer and components (or steps) of a system ofthe invention using a wireless network using a protocol such as Wi-Fi(IEEE standards 802.11, 802.11a, 802.11b, 802.11e, 802.11g, 802.11i, and802.11n, just to name a few examples). For example, signals from acomputer may be transferred, at least in part, wirelessly to componentsor other computers.

A specific type of electronic design automation tool is a circuitsimulation program or system. A circuit simulation program performsanalysis of circuits containing resistors, capacitors, inductors, mutualinductors, independent voltage and current sources, dependent sources,transmission lines, and semiconductor devices including diodes, bipolarjunction transistors (BJTs), junction field effect transistors (JFETs),and metal over semiconductor field effect transistors (MOSFETs). Acircuit simulator may perform nonlinear DC, nonlinear transient, linearAC, and other analyses.

One circuit simulation program is SPICE, originating from the Universityof California, Berkeley. SPICE stands for “Simulation Program IntegratedCircuits Especially!” Despite the success of SPICE and other circuitsimulation programs, existing circuit simulation programs usemethodology and computational techniques not suitable for use indistributed computing environment and determining real or exactsolutions.

FIG. 3 shows a black box model 304. In the figure, a block 304 has fourI/O (input/output) pins 308 a-308 d. The input-out pins of block 304 arethe points through which the block connects externally to other blocks(nodes and branches).

FIG. 4A shows a first Barycenter compact model block 404 with aninternal center point. The Barycenter compact model block has a boundarynode 408, a boundary branch 412, and a center point 416. The centerpoint can be one of the boundary points or internal node. The Barycenteris a coupling model between all boundary branches. The boundary branchdirection may be used +1 to represent current injection from boundarynode to the center point; −1 to represent current exit from center pointto boundary point. The Barycenter compact model may be a completecoupling modeling between all boundary nodes and boundary branches.

In this model block, there are a number of four boundary nodes 408 (orinput-output (I/O nodes). There are a four boundary branches (412) b1,b2, b3, and b4. Each of the four boundary nodes is connected to onebranch b1, b2, b3, or b4. There is a center point node 416 connected tobranches b1, b2, b3, and b4. Note that each the center point node isconnected to each of the boundary nodes through, at most, a singleboundary branch.

In comparison to other models for a circuit block, for the Barycentercompact model, all the boundary nodes connect through a single branch toa center point node. There are no boundary nodes that do not connect toother boundary nodes directly or connect to other nodes within theblock. There are no other nodes within the block other than the centerpoint node. A boundary node does not connect passing through anothernode before connecting to the center point node; the boundary nodedirectly connects to the center point node through a branch (e.g., b1,b2, b3, or b4).

Although referred to as the center point node, this node is not at aphysical center of the block. Center point is a reference that it is apoint between each of the boundary nodes. Center point 416 is internalto block 404, not one of the boundary or I/O nodes. The point is notnecessarily centered or in a center from a geometric point of view.

Each of the branches can include an impedance element such asresistance, capacitance, or inductance, or combination of these. Eachbranch can include a current or voltage source (e.g., independentcurrent source or independent voltage source) by itself or incombination with an impedance element.

Each branch impedance and source will have associated values orcoefficients. In a technique, a circuit (e.g., a circuit block or cell)is converted to a Barycenter compact model, which is an equivalent tothe circuit. There may be a number or family of equivalent models (notnecessarily Barycenter compact models). The circuit may have anyarbitrary interconnection of branches and nodes. However, after theconversion, the Barycenter compact model representation will have acircuit organization (which may be referred to as a “star pattern”) asshown or described in FIG. 4A (or FIG. 4B below). The impedance elementsand sources and coefficients for each of the elements and sources arenot the same as the original circuit, since the circuit graph isdifferent. During the conversion to the Barycenter compact model, thesecoefficients are calculated based on the original circuit.

FIG. 4B shows a second Barycenter compact model block 434 with aboundary node as a center point. Block 434 is an alternative to block404. For model 434, one of the boundary nodes or input-output (I/O) pinsis selected as the center point of the Barycenter compact model. In thespecific embodiment of FIG. 4B, boundary node 436 is selected as thecenter point. All remaining boundary nodes (438) of block 434 aredirectly connected to boundary node 436 using a boundary branch (432,e.g., b1, b2, and b3). The boundary branches connect to the boundarynode without passing through other boundary nodes or an internal centerpoint node.

Block 434 can be used to replace a block 404 configuation, and viceversa. For a hierarchical netlist, before simulation, the cells of thenetlist replaced with Barycenter compact model blocks, one model blockper cell, to improve the simulation speed. Each Barycenter compact modelblock is an equivalent model to the cell it is replacing. The top-levelof the netlist is also represented by a Barycenter compact model block,which incorporates lower-level Barycenter compact model blocks. TheBarycenter compact model for the hierarchical netlist may includeBarycenter compact model blocks with internal center point nodes,boundary node center points, and any combination of these. This allowsgreater flexibility in obtaining the Barycenter compact model for thenetlist, thus improving simulation speed and convergence.

Furthermore, due to their equivalence, a Barycenter compact model blockwith internal center point node can be replaced with a Barycentercompact model block with boundary node center points, and vice versa. Incertain simulation circumstances, one type of Barycenter compact modelblock may be favored over the other for performance reasons (e.g.,improves matrix convergence). Therefore, in a specific implementation, acell hierarchical circuit can be initially converted to a first type ofBarycenter compact model block (e.g., internal node center point) canlater be converted to a second type of Barycenter compact model block(e.g., boundary node center point), or vice versa. Such a change can bemade for certain circuit cells where it is recognized that the changewill improve performance.

FIG. 5A shows a black box model 504 being replaced or substituted with afirst Barycenter compact model 508 with an internal center point 509. Aswas discussed above, model 508 has boundary nodes, each of which isconnected through a single boundary branch to a single center pointnode. The single center point node is an internal node to the black box.

FIG. 5B shows a black box model being replaced with a second Barycentercompact model 511 with a boundary node as a center point. Model 511 isan alternative to model 508. For model 511, one of the boundary nodes orinput-output (I/O) pins is selected as the center point of theBarycenter compact model.

In this embodiment, a boundary node 513 is selected as the center point.A boundary node 521 is directly connected through a boundary branch b1to boundary node center point 513. A boundary node 524 is directlyconnected through a boundary branch b2 to boundary node center point513. A boundary node 526 is directly connected through a boundary branchb3 to boundary node center point 513. A boundary node 528 is directlyconnected through a boundary branch b4 to boundary node center point513.

In other Barycenter compact model embodiments for black box 504, anotherboundary node of the black box, such as boundary node 521, 524, 526, or528, can be selected as the center point instead of boundary node 513.Then the other boundary nodes will be directly connected to the selectedboundary node center point. The boundary nodes are the inputs, outputs,or input-outputs of the black box. The Barycenter compact models are anequivalent circuit to the black box and cell within the hierarchicalnetlist that the model is being used to replace.

In a technique of hierarchical circuit simulation, each black box modelis replaced with a Barycenter compact model, such as shown in FIGS. 5Aor 5B. Each Barycenter compact model is an equivalent circuit to thecell within the hierarchical netlist that the model is being usedreplacing. Since the netlist is hierarchical, within a black box orBarycenter compact model, there can be lower level black boxes orBarycenter compact models (just as there are lower level cells withinthe hierarchy).

A circuit simulation may include only Barycenter compact models withinternal center point nodes. A simulation may include only Barycentercompact models with boundary point center point nodes. A simulation mayinclude a combination of Barycenter compact models having internalcenter point nodes and boundary point center point nodes. In otherwords, within the same hierarchical circuit simulation, some of thecells are replaced with Barycenter compact models having internal centerpoint nodes and other cells are replaced with Barycenter compact modelshaving boundary pin center point nodes. The Barycenter compact modelsdiscussed are compatible with each other and can interoperate with eachother within the same circuit simulation. This gives greater flexibilitywhen performing circuit simulation.

FIG. 6 shows an implementation of interconnected blocks within a block.More specifically, a block A 604 is interconnected to a block B 608within a block C 610. Internal branches 612 a and 612 b connect blocks Aand B. A value calculated in block A can be transmitted to block B viainternal branch 612 a. Further, an output of block A can be received asinput to block B via internal branch 612 b. There is also a feedbackpath 614 from block B back to an I/O pin of block A. Block B has I/Opins 616 a and 616 b.

Block B is shown as having three I/O pins. However, a block can have anynumber of I/O pins. In other implementations, block B has one, two,five, or seven I/O pins. Further, in this specific implementation, twoblocks are interconnected. However, there can be more than twointerconnected blocks. For example, in another specific implementation,three, four, five, ten, or more blocks may be interconnected.

In an implementation, blocks A 604 and B 608 are substituted withBarycenter compact models, and the network is solved. FIG. 7 showsinterconnected blocks with Barycenter compact models for blocks A and B.A block A 704 (Barycenter compact model) and a block B 708 (Barycentercompact model) are interconnected within a block C 712. Internalbranches 716 a and 716 b connect blocks A and B at their center points.A value calculated in block A can be transmitted to block B via internalbranch 716 a. Further, an output of block A can be received as input toblock B via internal branch 716 b. There is also a feedback path 718from block B back to an I/O pin of block A. Connections between blocks Aand B may be referred to as intrablock branches because that arebranches or connections between the blocks.

Block B has I/O pins 720 a and 720 b. Block B is shown as having threeI/O pins. However, a block can have any number of I/O pins. In otherimplementations, block B has one, two, five, or seven I/O pins.

In an implementation, block A 604 is substituted for block A 704, andblock B 608 is substituted for block B 708.

FIG. 8 shows a block diagram of a Barycenter compact model with internalcenter point node. The Barycenter compact model compacts all internalnode information to boundary nodes and uses this model to represent thebox model. An advantage of this model is that it physically representsI/O boundary models of the black box. A block C 804 is generated withthe Barycenter compact model. Note that the internal branches andfeedback path can be modeled by (or collapsed into) a single Barycentercompact model. This facilitates ease in handling, making calculationssimpler, and allows for higher speed calculations.

Block C includes lower level blocks A and B. In this example, block C isa Barycenter compact model with internal center point node. However, asdiscussed above, block C could has be represented using a Barycentercompact model with a boundary node center point node and represented inFIGS. 4B and 4C.

FIG. 9 shows a graphical diagram of a system environment input. In animplementation, computers environment can be multicore, networked orgrid computing.

In a specific implementation, the database is a hierarchical database.An electronic circuit design may be provided as a hierarchical Spicenetlist. The user may provide source vectors for use in evaluating thedesign; and device models to give transistor, impedance, or otherphysical characteristics of components of the electronic design.

In an implementation, the user can supply a schematic capture, and thiscan be used to generate the circuit netlist. There is a layout file ordatabase, such as supplied as a GDS2 file. The layout provides thephysical geometries of the electronic circuit design. This layout isused to generate the masks for the manufacture of an integrated circuitwith the electronic circuit design. A layout viewer can view the layoutgraphically on a display. A layout extractor extracts the parasiticsfrom the layout and associates it with corresponding nodes in thenetlist. Other analysis tools may be used to view or evaluate theresults of an output of this system, such as using a three-dimensionallayout viewer.

With the netlist, this is converted by a convertor to database forhandling. Processing include hierarchical partitioning, input-output(I/O) slicing, and hierarchical scheduling. [82] FIG. 10 shows a flowfor a bottom up calculation. A specific implementation of a bottom upcalculation begins at a first step A 1004 and proceeds to a step 1008.At step 1008, a primitive cells job list is created. At a step 1012, ajob list is submitted and a Barycenter compact model is created. At astep 1016, blocks hierarchical dependence job list is created. At a step1020, a job list is submitted and a Barycenter compact model is created(same step as 1012). At a step 1024, the bottom calculation iscompleted. The output of these steps can be submitted to computersenvironment 1028.

FIG. 11 shows a flow for a top-down update. A specific implementation ofa top down update begins at a first step B 1104 and proceeds to step1108. At a step 1108, boundary branches, nodes, internal branches,nodes, and power of all blocks are updated. At a step 1112, boundarybranches, nodes, internal branches, nodes, and power of each primitiveare updated. At a step 1116, the top down update is completed. Theoutput of these steps can be submitted to computers environment 1120.

FIG. 12 shows a flow of an output. The flow starts at a first step 1204,and proceeds to a step 1208 in which data from a database 1208 isretrieved. In a specific implementation, the database is a hierarchicaldatabase. At a step 1212, a layout view extracts data from the database.At a step 1216, data is extracted from the database and a waveform isdisplayed by a waveform displayer. At a step 1220, data is extractedfrom the database and a schema is captured by a schematic capture. At astep 1224, data is extracted from the database and a layout is extractedby a layout extractor. Further, data from the database can betransmitted and received to and from computers environment 1228.

FIG. 13 shows a specific implementation of a hierarchical partition1301. The hierarchical partition has a top 1304, blocks 1308 a and 1308b, and primitive cells 1312 a and 1312 b. The hierarchical partition isshown as having four levels. However, the hierarchical partition canhave more than four levels or less than four levels. For example, thehierarchical partition can have three, five, seven, or eight levels.Further, a hierarchical partition can have any positive number of blocksor primitive cells. Any partitioning technique can be used to generate ahierarchical structure. Some examples include Binary, Quart Tree, K-dtree, Bin, and so forth.

FIG. 14 shows a hierarchical tree scheduler. The figure shows a specificimplementation of a hierarchical dependence tree 1401 with four levels.A top 1404 of the hierarchical dependence tree is at a level 4 (1408).As indicated, three blocks are at a level 3 (1408). As indicated, threeblocks 1412 are at a level 2 (1412). The remaining eleven primitivecells are at a level 1 (1416). There can be any number of primitivecells or blocks.

FIG. 15 shows a hierarchical dependence scheduler. At a level 4, block1504 transmits jobs to job list 1508 and can also transmit data to block1512. At a level 3, block 1512 transmits jobs to job list 1516 and canalso transmit data to block 1520. At a level 2, block 1520 transmitsjobs to job list 1524 and can also transmit data to block 1528. At alevel 1, block 1528 transmits jobs to job list 1532. In a specificimplementation, each job in the same level is fully independent.

In this scheduler, calculations occur in a top down fashion, so level 4first (by job list 1508), then level 3 (by job list 1516), then level 2(by job list 1524), and then level 1 (by job list 1532).

FIG. 16 shows a large I/O slicing into smaller I/O blocks. A block withlarger I/O 1604 slices a large I/O block into hierarchical less I/Oblocks. Block 1604 transmits data to a block A 1608, block B 1612, andblock C 1616, and block A 1608 can transmit data to block B 1612 whichcan transmit data to block C 1616. In an implementation, the usage ofmain memory size is controlled to fit the computer system memorylimitation.

FIG. 17 shows a flow diagram of a static simulation. The simulationbegins at a step 1704 in which a netlist database is provided. At a step1708, data from the netlist database is used as input. In animplementation, the input is a netlist parser, netlist hierarchicalpartition to generate primitive cells and blocks with an I/O boundary,netlist I/O slicing to create a next hierarchical level for a large I/Oboundary, or a hierarchical job scheduler and hierarchical data base. Ata step 1712, the input is used as input into a bottom up calculation.

During the bottom-up calculation, a primitive cells job list can becreated or primitive cells job list in parallel (multicores) ordistributed network computing can be submitted and the primitive cellsequivalent model using direct solver simultaneously can be generated.Further, a dependence job list of the blocks using a hierarchicalscheduler can be created, and a job list of blocks in parallel(multicores) or distributed network computing can be submitted and theequivalent I/O model using direct solver simultaneously can begenerated.

At a step 1720, a top down update is performed. In an implementation,during this step, boundary and internal branch current, nodal voltageand power of blocks in parallel (multicores) or distributed networkcomputing are updated simultaneously. Further, all boundary and internalbench current, bench voltage and power of primitive cells in parallel(multicores) or distributed network computing are updatedsimultaneously. In other implementations, they are not updatedsimultaneously.

At a step 1724, a result of the top down update is stored in a database.At a step 1728, data is extracted from the database and is output.Output can be branch currents, branch voltages, nodal voltages, blockpower, or branch power. Further, output can be a waveform display,layout mapping display, highlight EM problems, highlight static anddynamic IR drop, signal noise, signal timing, signals cross talk, andpower consuming. Computers environment 1732 can be multicores andnetworked.

FIG. 18 shows a flow diagram of a dynamic simulator.

FIG. 19 shows a flow diagram of a circuit simulator. The simulationbegins at a step 1904 in which a netlist database is provided. At a step1908, data from the netlist database is used as input. At a step 1912,vectors from a database 1916, device models from a database 1920, andupdates from a top down update are used to perform a bottom upcalculation. At a step 1924, an output of the bottom up calculation isused to perform a top down update. At a step 1928, a result from the topdown update is stored in a database. At a step 1932, data is extractedfrom the database and is output.

Using the Barycenter compact model and hierarchical scheduler, thenonlinear device model can be integrated into the system. Further, usingthe Newton-Raphson method, the nonlinear boundary problems can besolved.

In an implementation, the invention has several advantages, such as nomemory sharing, ability to control the usage of memory size, fullyutilizing all CPUs power linearly, ease in programming in a parallel anddistributed computing system environment, faster performance andcapacity, and if a direct matrix solver is used, the exact answer can becalculated regardless of the different partitions.

U.S. patent application Ser. Nos. 11/421,206 and 11/421,212, both filedMay 31, 2006, describe network tearing and global and local links. Thenetwork tearing techniques described in those patent applications can beused in conjunction with the techniques described in this patentapplication to solve circuit networks of hierarchical electronic circuitdesigns.

FIG. 20 shows an example of a branch of a circuit. This branch has acurrent source 2002, impedance or admittance 2005, and voltage source2009. In electrical engineering, the admittance (Y) is the inverse orreciprocal of impedance (Z). There are also currents I_(b) and J_(b).

A circuit component is represented using a branch such as shown in FIG.7. More specifically, each tree branch representsresistance-capacitance-inductance (RCL) and one or more sources, such ascurrent source, voltage source, dependent current source, and dependentvoltage source. As discussed above, a device model for a transistor orother device is a RCL network with a dependent source. Therefore, eachtransistor of the circuit will be a branch in the network graph. Theentire graph is a RCL network with sources.

FIG. 21 shows a graph of a tree and links. Tree branches are shown usingsolid lines, and links are shown using dotted lines. At a junction (orintersection) of two or more tree branches is a node. In a graph, someinformation is associated with each node and edge. For example, a nodeof the graph may be a node in the circuit and edge may be a branch inthe circuit. A graph is an abstract data type that consists of a set ofnodes and a set of edges that establish relationships or connectionsbetween the nodes. A specific technique, among others, to implement agraph data structure is to use linked lists over the nodes and theingoing and outgoing edges of the nodes. Another technique to implementa graph is to use an array structure. There are many other approaches toimplementing a graph data structure in a computer system, and any ofthese may be used.

After a graph of the circuit is built, a technique identifies branchesand links of the tree. A tree of a graph is a connected subgraph thatincludes all the nodes of the graph but contains no loops. A loop is acircular path from a first node through other nodes and returns to thefirst node. A loop has a voltage drop of zero. A subgraph of a graph isa set of branches and nodes belonging to a graph. A link forms a loopwith one or more tree branches, and may be referred to as a unit link. Atree branch does not form a loop. Typically, a circuit or system has onetree and multiple links. FIG. 21 shows an example of a graph with nodes,branches, and links A circuit network is represented using such a graphdata structure.

For a hierarchical circuit design, a node of the graph may be a subcellor subcircuit of the hierarchical circuit design. Nodes at the lowestlevel (or bottom level) of the hierarchy may be referred to asprimitives or leafs.

There are many techniques to identify trees and links One approach is adepth first search. Another technique is a breadth first search. Eitherof these may be used or a combination of these two may be used. Forexample, one technique of finding a tree involves starting at a startingnode. This node can be any node in the graph. Depending on which node isselected as the starting node, the tree may be different from a treefound using a different starting node. Proceeding branch by branchthrough the graph, each branch of the graph will be designated as a treeor link to create a graph, such as in FIG. 8, having tree branches andlinks.

Partitioning the tree into subtrees breaks up or tears a tree intosmaller subtrees. In an embodiment of the invention, partitioning thetree breaks up the tree into a number of subtrees with the same or closeto the same number of tree branches. Each subtree is a subcircuit of thecomplete circuit. This technique may be automated using, for example, acomputer. In particular, a number count or predefined count may beselected, such as two, three, four, five, six, seven, eight, ten, morethan ten, and so forth. Smaller counts are used to break the tree intomore subtrees than for larger counts.

According to one technique, a first subtree is formed by repeatedlyadding branches of the tree to the first subtree when a number ofbranches in the first subtree is less than a predefined count and thereare branches in the tree which have not yet been assigned to a subtree.This technique is continued with the second subtree, third subtree, andso forth until all branches of the tree have been assigned to a subtree.Each branch of the tree is only assigned to one subtree. Each subtree isa tree with connected branches. Using this counting approach, eachsubtree will have a number of branches less than or equal to thepredefined count. The predefined count may be set by the system, or maybe user defined. This technique can be used for flat electronic designsor designs where hierarchy need not be maintained (where the hierarchycan be flattened).

For a hierarchical circuit design, an approach to partitioning the treeinto subtrees (while maintaining the hierarchy) is to recursivelypartition the tree into smaller and smaller subtrees while a number ofbranches of each succeeding subtree is greater than a predefinedpartition size. More specifically, a tree is partitioned into a numberof subtrees. The size of each subtree is checked against the predefinedpartition size. If the size of a subtree is larger than the predefinedpartition size (e.g., greater number of branches than the predefinedpartition size), that subtree is partitioned again. This technique isapplied recursively for each succeeding subtree obtained, until each ofthe resulting subtrees is equal to or smaller than the predefinedpartition size.

The predefined partition size may be defined by the user, such as in aparameter file read by software before or during runtime. Recursivepartitioning maintains the hierarchy in a system graph of the electroniccircuit design. In an implementation, the recursive partitioning routinecontinues until the subtree size is smaller than the defined partitionsize or close to the predefined partition size in case when a primitivecell is obtained that cannot be further partitioned.

In another implementation, an initial system graph is divided or brokeninto a number of instances before the recursive partitioning. Eachinstance is a subtree. The number of instances is predefined, and can beuser defined similarly to the predefined partition size. The number ofinstances may be related to the number of processors or processing coresavailable. For example, for a quad-core processor, the number ofinstances selected may be four. For eight cores (available on a singlemachine or multiple machines together), the number of instances may be8. For a 64-core system, the number of instances may be 64, and soforth. Then, each of the instances (which is a subtree) is thenrecursively partitioned while a number of branches of each succeedingsubtree is greater than a predefined partition size.

For a hierarchical circuit design, nodes at the lowest level (or bottomlevel) of the hierarchy may be referred to as primitives or leafs. In atypical hierarchical circuit design, there can be any number ofprimitive cells, and there are multiple instances of a primitive cell.

For example, FIG. 22 shows the tree of FIG. 21 partitioned into threesubtrees (also may be called local trees or subblocks), where thesubtrees have 8 or 9 branches. As a further example, if the count istwo, a tree is separated every two tree branches until it can no longerbe subdivided. If there are an odd number of branches in the tree, thenthe last subtree will have one branch, rather than two.

Described above is merely an example of some techniques of partitioning.Other techniques may be used to partition the tree into subtrees, andany of these other techniques may also be used in implementing theinvention.

A technique determines which links are local are which are global. InFIG. 22, tree branches are shown using solid lines, global links areshown using broken lines, and local links are shown using dotted lines.A link that forms a loop in a subtree is a local link. A link that formsa loop in two or more subtrees (or multiple subtrees) is a global link.A global link forms a loop with branches of two or more subtrees.

FIG. 9 shows a tree having three subblocks. For ease in identifyingsubblocks in the figure, each subblock is circled with a boundary line.Each subblock has a local tree and local links. The local links for aparticular subblock are located within the subblock and do not cross asubblock boundary line. There are global links between subblocks andlocal trees. The global links cross the subblock boundary lines.

A subblock solver will solve each subblock including local links, butnot global links. During this step, a circuit subblock is simulatedwhile circuit branches between subblocks (i.e., global links) areignored. The subblock solver will take individual subtrees or subblocksand solves each of these subtrees or subblocks independently.

Therefore, a solution to any particular subblock may be determinedwithout considering other subblocks. So, it will be immaterial in whatthe order the subblocks are simulated. This helps a distributedcomputing embodiment of the invention because by allowing each subblockto be simulated independently, there will not be any schedulingproblems.

The solver determines the IR drop (current-resistance drop) or voltagedrop for each branch or local link, or both. In other words, the solverdetermines the voltage and current for each branch for the givenconditions. To determine the IR drops, the computer will perform amatrix inversion or LU (lower-upper) factoring of a subblock in the treenetwork. The subblock which is solved includes branches and local links,but not global links.

With an embodiment of a system of the invention, a single computer mayperform computations for each subblocks in sequence. Since the subblockis much smaller than the matrix of the entire tree, computation time isreduced. In further embodiments, multiple computers may performcomputations for subblocks at the same time, reducing simulation timecompared to using one computer or solving a very large matrix for thewhole tree.

More specifically, the subblock solver can send each of the subtrees asa computing task or job to be solved by a different computer in adistributed computing network. If there are too few computers to do alljobs simultaneously, two or more tasks may be queued up on a computerfor computer. In such fashion, different computers in the distributedcomputing network can perform parts of circuit simulation calculationsin parallel, thus further speeding up the circuit simulation process.

An interblock solver will solve the IR drop or voltage drop for theglobal links, which are the links or interconnection between the blocks.The interblock solver determines the junction voltage or junctioncurrent, or both, for each global link. The junction voltage is thevoltage across the global link, and the junction current is the currentwhich flows through the global link. The results of the interblocksolver will be independent of the results of subblock solver. So, theinterblock solver does not use the results from the subblock solver.Interblock solver may be performed using a single computer ordistributed computing.

A technique combines the partial results with the interconnect-levelresults (obtained using the interblock solver) to find the exact or realresults for the entire circuit. Update system partial results updatesthe system results with the contribution of the global links to eachsubtree. In an implementation, the results obtained after update systempartial results will be the real solution for the whole system, as ifthe system were solved together as one large matrix. This real solutionwill not be an estimation, approximation, or an iteratively obtainedsolution, but an exact solution.

An output block outputs the IR drop results into database. For eachnode, there may be a nodal voltage, branch voltage, and branch current.Branch voltage is the voltage across the two nodes of a branch. Giventhe nodal voltages, the branch voltage may be calculated. Branch currentis the current flowing through a branch.

Some specific flows for circuit simulation are presented in this patent,but it should be understood that the invention is not limited to thespecific flow and steps presented. A flow of the invention may haveadditional steps (not necessarily described in this application),different steps which replace some of the steps presented, fewer stepsor a subset of the steps presented, or steps in a different order thanpresented, or any combination of these. Further, the steps in otherimplementations of the invention may not be exactly the same as thesteps presented and may be modified or altered as appropriate for aparticular application or based on the data

In an implementation, a system includes:

1. A network.

2. A number of computing devices, connected to the network.

3. A circuit simulation block including:

3a. A system tearing block to identify each branch of a givenresistance-capacitance-inductance (RCL) circuit network graph as atleast one of a tree branch or a link, where the tree branches form atree of the graph, and the tree includes no loops.

3b. Divide the tree into at least a first subtree instance and a secondsubtree instance, where each the first and second subtree instancesincludes no loops, and

3c. Identify links of the given RCL circuit network graph as at leastone of a global link or a local link, where a global link forms a pathfrom a branch of the first subtree to a branch of the second subtree, aglobal link forms a loop with branches of the first and second subtrees,the first subtree is recursively partitioned to obtain i subtrees whilea number of branches of each of the i subtrees is greater than apredefined partition size, and the loop with the global link andbranches of the first and second subtrees has a voltage drop of 0; and

3d. A subblock solver block to send an i subtree, without any globallinks, through the network to a first computing device of the pluralityof computing devices for calculation and an i+1 subtree, without anyglobal links, through the network to a second computing device of theplurality of computing devices for calculation.

Further, the second subtree can be recursively partitioned to obtain jsubtrees while a number of branches of each of the j subtrees is greaterthan the predefined partition size. The subblock solver block is to senda j subtree, without any global links, through the network to a thirdcomputing device of the computing devices for calculation and a j+1subtree, without any global links, through the network to a fourthcomputing device of the plurality of computing devices for calculation.

Each of the i subtrees is modeled by a block including: boundary nodes;boundary branches, one boundary branch connected to each boundary node;and a center point node connected to each of the of the boundary nodesthrough, at most, a single boundary branch.

In an implementation, a method includes:

1. Recursively partitioning a system graph of an hierarchical electroniccircuit design into subcells while a number of branches in each of thesubcell is greater than a predefined partition size, where subcells at alowest level in a tree of the system graph are referred to as primitivecells. The hierarchical electronic circuit design may be specified usinga hierarchical netlist.

2. Converting each of the subcells into corresponding model blocks, eachmodel block including: boundary nodes; boundary branches, one boundarybranch connected to each boundary node; and a center point nodeconnected to each of the of the boundary nodes through, at most, asingle boundary branch.

3. Determining initial value solutions for primitive cells using modelblocks for the primitive cells.

4. Assembling the primitive cells in a hierarchical fashion from abottom to a top of the tree to obtain an assembled structure whichmaintain the hierarchy of the hierarchical electronic circuit design.

5. Determining boundary conditions for the subcells.

6. Using the determined initial value solutions for the primitive cells,calculating initial value solutions for each subcell from a bottom to atop of the assembled structure.

8. Using the determined boundary conditions, calculating boundary valuesolutions for each subcell from a top to a bottom of the assembledstructure.

The calculating initial value solutions for each subcell from a bottomto a top of the assembled structure can include: calculating a solutionto each subcell while ignoring global links, where a global link forms apath from a branch of a first subcell to a branch of a second subcell,and a loop with the global link and branches of the first and secondsubcells has a voltage drop of 0.

The method may further include building an input-output connectivity ofthe assembled structure.

The assembling the primitive cells in a hierarchical fashion from abottom to a top of the tree to obtain an assembled structure mayinclude:

1. Providing a first model block, including first boundary nodesconnected via first boundary branches to a first center point node,corresponding to a first subcell.

2. Providing a second model block, including second boundary nodesconnected via second boundary branches to a second center point node,corresponding to a second subcell;

3. Determining a hierarchical connectivity between the first and secondmodel blocks from the system graph tree includes an intrablock branch,the intrablock branch including at least one of:

3a. A first intrablock branch connecting the first center point node ofthe first model block to the second center point node of the secondcenter point node,

3b. A second intrablock branch connecting one of the first boundarybranches of the first model block to one of the second boundary branchesof the second model block, OR

3c. A third intrablock branch connecting the second center point node ofthe second model block to one of the first boundary nodes of the firstmodel block.

4. Obtaining a third model block corresponding to the first and secondmodel blocks and determined hierarchical connectivity, where the thirdmodel block includes third boundary nodes connected via third boundarybranches to a third center point node, where the third center point nodeis connected to each of the third boundary nodes through, at most, asingle third boundary branch.

In an implementation, a method includes:

1. Providing a hierarchical system graph of an electronic circuitincluding multiple instances of at least one primitive.

2. Identifying each branch of the system graph as at least one of a treebranch or a link, where the tree branches form a tree of the systemgraph without flattening a hierarchy of the system graph.

3. Dividing the tree into n subtree instances including at least firstand second subtree instances, where n is predefined number of instanceshaving an integer value 2 or greater.

4. Identifying links of the system graph as at least one of a globallink or a local link, where the system graph includes both global andlocal links, a global link forms a path from a branch of the firstsubtree instance to a branch of the second subtree instance, and a loopwith the global link and branches of the first and second subtrees has avoltage drop of 0;

5. Recursively partitioning the first subtree instance into i subtreeswhile a number of branches of each of the i subtrees is greater than apredefined partition size, where the partition size is definedseparately from the number of instances.

6. Identifying links forming a path between branches of the i subtreesof the first substree instance as global links.

7. Recursively partitioning the second subtree instance into j subtreeswhile a number of branches of each of the j subtrees is greater than thepredefined partition size.

8. Identifying links forming a path between branches of the j subtreesof the second substree instance as global links.

The recursive partitioning the first subtree instance into i subtreeswhile a number of branches of each of the i subtrees is greater than apredefined partition size does not include flattening the hierarchy ofthe system graph. Therefore, the hierarchy of the system graph of theelectronic circuit is maintained. This speeds up calculations becauseprimitive cells need only be simulated once and then assembled together.When flattened, the hierarchy is lost and instances of the primitivesare recalculated each time. Some primitives, such as an inverter, may berepeated many times in an electronic or integrated circuit design.Therefore, in an implementation, each of the i subtrees of the firstsubtree instance and j subtrees of the second subtree instance maintainsthe hierarchy of the system graph of the electronic circuit.

Further, calculating a solution to each of the i subtrees initiallyignores global links, and calculating a solution to each of the jsubtrees initially ignores global links. Calculating a solution to oneof the i subtrees is performed on a different computing device thancalculating a solution to one of the j subtrees. This facilitatesdistributed computing. Calculating a solution to each of the i subtreeswhile ignoring global links can include performing a matrix inversion.Calculating a solution to each of the i subtrees while ignoring globallinks can include performing LU factoring.

A local link forms a loop in one of the i subtrees of the first subtreeinstance. A tree branch forms no loops in the system graph. A link formsa loop in the system graph. Each of the i subtrees is modeled by a blockincluding: boundary nodes; boundary branches, one boundary branchconnected to each boundary node; and a center point node connected toeach of the boundary nodes through, at most, a single boundary branch.

The system may optionally include a self-verification tool. Theself-verification tool may be included in some embodiments, while it isomitted in other embodiments of the invention. The self-verificationtakes the simulation results of verifies the results are correct. U.S.patent application Ser. No. 11/279,391, filed Apr. 11, 2006, issued asU.S. Pat. No. 7,461,360 on Dec. 2, 2008, discusses simulationverification and is incorporated by reference along with any otherreferences cited in this application. A technique of simulationverification involves determining whether the simulation resultsobtained satisfy Kirchhoff's current law (KCL), Kirchhoff's voltage law(KVL), and the power conservation law. If any of these three laws arenot satisfied, a not-verified condition results, which means thecalculated simulation results contain an erroneous. For example, if KCLis verified, but not KVL and not power, this results in a not-verifiedcondition. If KCL and KVL are verified, but not power, this results in anot-verified condition. In an implementation, to obtain a verifiedcondition, KCL, KVL, power must be verified, otherwise a not-verifiedcondition results.

For example, a method may further include:

1. Using the calculated initial value and boundary value solutions,determining node voltages for each node of the hierarchical electroniccircuit design.

2. Using the node voltages, determining branch voltages for branches inthe system graph.

3. Identifying independent loops in the graph.

4. Summing the voltages for each independent loop in the graph.

5. Summing the currents at each node in the graph.

6. Summing power consumed for each branch in the graph to obtain a totalpower consumed.

7. Determining a total input power to the circuit network using theinput sources associated with the circuit.

8. Subtracting the total power consumed from the total input power toobtain a total power difference.

9. Indicating a not-validated condition when at least one of the loopsin the graph has a nonzero sum;

10. Indicating a not-validated condition when at least one of the nodesin the graph has a nonzero sum.

11. Indicating a not-validated condition when the total power differenceis not zero.

In another implementation, a technique includes: providing a circuitnetwork specified in a netlist format and input sources associated withthe circuit; providing a simulation output for the circuit, where thesimulation output includes node voltages for each node of the circuit;building a graph data structure from the circuit netlist; using the nodevoltages, determining branch voltages for branches in the graph;identifying a tree and links in the graph; identifying independent loopsin the graph; summing the voltages for each independent loop in thegraph; summing the currents at each node in the graph; summing powerconsumed for each branch in the graph to obtain a total power consumed;determining a total input power to the circuit network using the inputsources associated with the circuit; subtracting the total powerconsumed from the total input power to obtain a total power difference;indicating a not validated condition when at least one of the loops inthe graph has a nonzero sum; indicating a not validated condition whenat least one of the nodes in the graph has a nonzero sum; and indicatinga not validated condition when the total power difference is not zero.

In various implementations, a current for a branch is calculated bybranch voltage divided by impedance for the branch. A validatedcondition is indicated when every loop in the graph sums to zero, everynode in the graph sums to zero, and the total power difference is zero.In a graphical viewer, nodes of the circuit are highlighted wherenonzero current summing results were obtained. In a graphical viewer,branches of the circuit are highlighted where nonzero voltage summingresults were obtained.

This description of the invention has been presented for the purposes ofillustration and description. It is not intended to be exhaustive or tolimit the invention to the precise form described, and manymodifications and variations are possible in light of the teachingabove. The embodiments were chosen and described in order to bestexplain the principles of the invention and its practical applications.This description will enable others skilled in the art to best utilizeand practice the invention in various embodiments and with variousmodifications as are suited to a particular use. The scope of theinvention is defined by the following claims.

The invention claimed is:
 1. A method comprising: providing ahierarchical system graph of an electronic circuit comprising multipleinstances of at least one primitive; identifying each branch of thesystem graph as at least one of a tree branch or a link, wherein thetree branches form a tree of the system graph without flattening ahierarchy of the system graph; dividing the tree into n subtreeinstances comprising at least first and second subtree instances,wherein n is predefined number of instances having an integer value 2 orgreater; identifying links of the system graph as at least one of aglobal link or a local link, wherein the system graph comprises bothglobal and local links, a global link forms a path from a branch of thefirst subtree instance to a branch of the second subtree instance, and aloop with the global link and branches of the first and second subtreeshas a voltage drop of 0; using at least one computer processor,recursively partitioning the first subtree instance into i subtreeswhile a number of branches of each of the i subtrees is greater than apredefined partition size; modeling each of the i subtrees using a firstblock model or a second block model, wherein the first block modelcomprises: a plurality of boundary nodes; a plurality of boundarybranches, one boundary branch coupled to each boundary node; and aninternal point node coupled to each of the boundary nodes through, atmost, a single boundary branch
 2. The method of claim 1 wherein thepartition size is defined separately from the number of instances. 3.The method of claim 1 wherein the identifying each branch of the systemgraph as at least one of a tree branch or a link, wherein the treebranches form a tree of the system graph comprising: not flattening ahierarchy of the system graph.
 4. The method of claim 1 comprising:identifying links forming a path between branches of the i subtrees ofthe first substree instance as global links; using at least one computerprocessor, recursively partitioning the second subtree instance into jsubtrees while a number of branches of each of the j subtrees is greaterthan the predefined partition size; and identifying links forming a pathbetween branches of the j subtrees of the second substree instance asglobal links.
 5. The method of claim 1 wherein the second block modelcomprises: an internal point node, which is a boundary node of thesecond block model; a plurality of boundary nodes which are not theinternal point node of the second block model; and a plurality ofboundary branches, one boundary branch coupling to each boundary node,which is not the internal point node, to the internal point node.
 6. Themethod of claim 1 wherein the recursively partitioning the first subtreeinstance into i subtrees while a number of branches of each of the isubtrees is greater than a predefined partition size comprises: notflattening the hierarchy of the system graph.
 7. The method of claim 2wherein each of the i subtrees of the first subtree instance and jsubtrees of the second subtree instance maintains the hierarchy of thesystem graph of the electronic circuit.
 8. The method of claim 2 furthercomprising: calculating a solution to each of the i subtrees whileignoring global links; and calculating a solution to each of the jsubtrees while ignoring global links.
 9. The method of claim 8 whereincalculating a solution to one of the i subtrees is performed on adifferent computing device than calculating a solution to one of the jsubtrees.
 10. The method of claim 8 wherein the calculating a solutionto each of the i subtrees while ignoring global links comprisesperforming a matrix inversion.
 11. The method of claim 8 wherein thecalculating a solution to each of the i subtrees while ignoring globallinks comprises performing lower-upper (LU) factoring.
 12. The method ofclaim 1 wherein a local link forms a loop in one of the i subtrees ofthe first subtree instance.
 13. The method of claim 1 wherein a treebranch forms no loops in the system graph.
 14. The method of claim 1wherein a link forms a loop in the system graph.
 15. The method of claim1 wherein a second block model incorporates a plurality of lower-levelfirst model blocks.
 16. A system comprising: a network; a plurality ofcomputing devices, coupled to the network; and a circuit simulationcomponent comprising: a system tearing component to identify each branchof a given resistance-capacitance-inductance (RCL) circuit network graphas at least one of a tree branch or a link, wherein the tree branchesform a tree of the graph, and the tree comprises no loops, divide thetree into at least a first subtree instance and a second subtreeinstance, wherein each the first and second subtree instances comprisesno loops, and identify links of the given RCL circuit network graph asat least one of a global link or a local link, wherein a global linkforms a path from a branch of the first subtree to a branch of thesecond subtree, a global link forms a loop with branches of the firstand second subtrees, the first subtree is recursively partitioned toobtain i subtrees while a number of branches of each of the i subtreesis greater than a predefined partition size, and the loop with theglobal link and branches of the first and second subtrees has a voltagedrop of 0; and a model replacement component to replace cells of thegiven RCL circuit network graph with a model, each model having at mosta single internal point node to which boundary nodes of the block arecoupled to through at most a single boundary branch.
 17. A methodcomprising: using at least one computer processor, recursivelypartitioning a system graph of an hierarchical electronic circuit designinto a plurality of subcells while a number of branches in each of thesubcell is greater than a predefined partition size, wherein subcells ata lowest level in a tree of the system graph are referred to asprimitive cells; converting each of the subcells into correspondingmodel blocks, each model block having only one node to which eachboundary node of the model block is directly coupled to; determininginitial value solutions for primitive cells using model blocks for theprimitive cells; determining boundary conditions for the subcells; andusing the determined initial value solutions for the primitive cells,calculating initial value solutions for each subcell from a bottom to atop of the assembled structure.
 18. The method of claim 17 comprising:assembling the primitive cells in a hierarchical fashion from a bottomto a top of the tree to obtain an assembled structure which maintain thehierarchy of the hierarchical electronic circuit design.
 19. The methodof claim 17 comprising: using the determined boundary conditions,calculating boundary value solutions for each subcell from a top to abottom of the assembled structure.
 20. The method of claim 17comprising: using the calculated initial value and boundary valuesolutions, determining node voltages for each node of the electroniccircuit design; using the node voltages, determining branch voltages forbranches in the system graph; and identifying independent loops in thegraph.